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A modified version of the conditional symmetry method, together with the 
classical method, is used to obtain new classes of elliptic solutions of the isentropic 
ideal compressible fluid flow in (3+1) dimensions. We focus on those types of 
solutions which are expressed in terms of the Weierstrass p-functions of Riemann 
invariants. These solutions are of special interest since we show that they remain 
bounded even when these invariants admit the gradient catastrophe. We describe 
in detail a procedure for constructing such classes of solutions. Finally, we present 
several examples of an application of our approach which includes bumps, kinks 
and multi-wave solutions. 

1 Introduction 

The purpose of this paper is to construct bounded elliptic solutions of a compressible 
isentropic ideal flow in (3 + 1) dimensions. Such solutions exist even in the case where 
the Riemann invariants admit the gradient catastrophe. 

Let us flrst present a brief description of a procedure detailed in [lOj for constructing 
rank-/c solutions in terms of Riemann invariants for the case of an isentropic compress- 
ible ideal fluid in (3 + 1) dimensions. Such a model is governed by the equations 
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a = 1,2,3,4, 



(1.1) 
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where A^.A"^ and A^ are 4x4 real- valued matrix functions of the form 
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J = 1,2,3, 



K = 2(7 — 1)^^ and 7 is the adiabatic exponent of the medium under consideration. 
The independent and dependent variables are denoted by x = (t = , , x'^ , x^) G 
X C M"^ and u = {a,u) E U G M^, respectively, and Ut stands for the first order partial 
derivatives of u, i.e. = du"'/dx\ a = 1,...,4, i = 0,1,2,3. Here, the quantity a 
stands for the velocity of sound in the medium and u is the velocity vector field of the 
flow. Throughout this paper, we adopt the summation convention over repeated lower 
and upper indices. 



The purpose of this article is to obtain rank-/c solutions of system (1.1) expressible 
in terms of Riemann invariants. To this end, we seek solutions u{x) of (1.1) defined 
implicitly by the following set of relations between the variables u", and x*, 

u = f{r\x,u),...,r\x,u)), r^ix,u) = Xfiu)x\ kei {X^I^ + A' (u) Xf) ^ , (1.2) 

for some function / : M'^ ^ and A = 1, . . . ,k < 3. Such a solution is called 
a rank-fc solution if rank(u°) = k. The functions r^{x,u) are called the Rie man n 
invariants associated with the wave vectors = (A;^, A^) G of the system ( |l.l[ ). 
Here, A"^ = (A^, A^, A3 ) denotes a direction of wave propagation and the eigenvalue Aq 
is a phase velocity of the considered wave. Two types of admissible wave vectors for 
the isentropic equations (1.1) are obtained by solving the dispersion relation 

det (Ao(m)J4 + Ai(M)^^(M)) = [{Xo + u ■ X)^ - a^X^]{Xo + u ■ Xf = 0. (1.3) 



They are called the entropic (E) and acoustic (S) wave vectors and are defined by 

i) X^ = {ea + u ■ e, —e), e = ±1, ii) A'^ = (det (m, e, m), — e x m), |e|^ = l, (1.4) 

where e and m are unit and arbitrary vectors, respectively. 

The construction of rank-/c solutions through the conditional symmetry method 
(CSM) is achieved by considering an overdetermined system, consisting of the original 



system (1.1) in 4 independent variables together with a set of compatible first order 



differential constraints (DCs), 

C(«)< = 0, AfC = 0, 



A-k, 



;i.5) 



for which a symmetry criterion is automatically satisfied. Such notions as conditional 
symmetry, conditional symmetry algebra and conditionally invariant solution for the 



original system (1.1) we use in accordance with definitions given in [T^ . Under the 
above circumstances, the following result holds : 
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The isentropic compressible ideal fluid equations (1.1) admit a (4 — /c)-dimensional 
conditional symmetry algebra L if and only if there exists a set of (4 — k) linearly 
independent vector fields 



Xa = Cain) 



d_ 



1,...,4-A;, ker(^'(n)Af) ^0, Af C = 0, A = 1, . . . , A; < 3, 

(1.6) 



which satisfy on some neighborhood of (xq, uq) E X x U the trace conditions 

= 0, fi = l, 



i) trj^^S^A 



dr 



0, ii) tr(^-^^^ ^ 



dr ^^""^ dr 



,4, 



;i-7) 



where 



A = (Af ) G I 
'dX^ 



du"^ 



pfcx4 



(r\...,r'=) eM^ 
s = 1, . . . , A; — 1, 



9r 



n4xfc 



and (ai, . . . , as) denotes the symmetrization over all indices in the bracket. Solutions of 
the system which are invariant under the Lie algebra L are precisely rank-fc solutions 
of the form (1.2). 

This result is a special case of the proposition in [TU]. Note that these symmetries 
are not symmetries of the original system, but they can be used to construct solutions 
of the overdetermined system composed of (1.1) and (1.5). 

For the case of rank-1 entropic solution E, the wave vector A^ is a non-zero multiple 
of ( 1.4 i). Therefore, the corresponding vector fields Xj and Riemann invariant r become 



Xi = —(a + e-u) 



1 ^ A 



1,2,3, 



r{x, u) = {a + u ■ e)t — e ■ X, \e\' 



where we chose e = 1 in (1.4 i). Rank-1 solutions invariant under the vector fields 
{Xi, X2, X^} are obtained through the change of coordinates 



t = t, x^ = r{x, u), x"^ = x'^, x^ = x^, a = a, v} = v} , = u^, = u'^, 



1.9) 



on X M^. Assuming that the direction of the wave vector e is constant, the fluid 
dynamics equations (1.1) transform into the system 

du^ 



da 



du^ 



dx^ * dx^ ' 

with the invariance conditions 



dx 



da 
^ dx^ ' 



1,2,3, 



;i.io) 



0, uf = 



x3 



0, J = 2,3, « = 1,2,3. 



(1.11) 
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The general rank-1 entropic E solution takes the form 



a{t, x) = a{x^), u\t, x) = Keia{x^) + Ci, Cj G M, « = 1,2,3, 



(1.12) 



where the Riemann invariant x^ = r{x,u) is given by 

r{x, u) = [(1 + K)a + e- C]t - e ■ x, C= {Ci, C2, C3) G Ml 
A similar procedure can be applied to the rank-1 acoustic solution S. Here, the 



wave vector A"^ is a non-zero multiple of (1.4 ii) and the corresponding vector fields Xj 
and Riemann invariant are 



r[x,u) 



det (u,e,m)t — (ex fh) ■ X, Xj 



(e X m)i d d 
det (u, e, rh) dt dx^ 



1,2,3. (1.13) 



Again, the change of variables (1.9) leads to transformed dynamical equations, which 



we integrate in order to find rank-1 acoustic solution of the form 

a(t, x)=ao, u^(t,x) = u^{x^), u^(t,x) = u^{x^), C G M, 

u^{t, x) = {eim2 - 62^1)"^ [C - {621713 - e3m2)u^{x^) - {esuii - eim^)^'^ (x^ 



(1.14) 



Here v} et v? are arbitrary functions of the Riemann invariant x^ = r{x,u) which has 
the explicit form 



r[x, u 



Ct — det (x, e, fh) 



In general, the overdetermined system composed of (1.7 i) and (1.7 



(1.15) 
ii) is nonlinear 

and cannot always be solved in a closed form. Nevertheless, particular rank-/c solutions 
for many physically interesting systems of PDEs are well worth pursuing. These par- 
ticular solutions of (1.7 i) and (1.7 ii) can be obtained by assuming that the function 



/ is in the form of a rational function, which may also be interpreted as a truncated 
Laurent series in the variables r^. This method can work only for equations having the 
Painleve property Consequently, these equations can be very often integrated in 
terms of known functions. 

Applying a version of the conditional symmetry method to the isentropic model 



(1.1), several new classes of solutions have been constructed in a closed form P ITO] . 
Comparing these results with the ones obtained via the generalized method of carac- 
teristics (CMC) [20j, it was shown that more diverse classes of solutions are involved 
in superpositions (i.e. rank-fc solutions) than in the case of the CMC [TU] . 

This paper is a continuation of the papers [HI ID]- The objective is to construct 



bounded elliptic solutions of the isentropic system (1.1) using the version of the CSM 



proposed in [lOj. These types of solutions are obtained through a proper selection 



of differential constraints (DCs) compatible with the initial system of equations (1.1). 



That is, the solution should satisfy both the initial system (1.1) and the differential 



constraints (1.5). Among the new results obtained, we have rank-2 and rank-3 peri- 



odic bounded solutions expressed in terms of Weierstrass p-functions. They represent 
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bumps, kinks and multiple waves, all of which depend on Riemann invariants. These 
solutions remain bounded even when the invariants admit a gradient catastrophe. 

The paper is organized as follows. In Section 2 we construct rank-2 and rank- 
3 elliptic solutions of the system, among which multiple waves and doubly periodic 
solutions are included, and we show that they remain bounded everywhere. Section 3 
summarizes the results obtained and contains some suggestions for future developments. 



2 Rank-2 and rank-3 solutions 

The construction approach outlined in Section 1 has been applied to the isentropic 



flow equations (1.1) in order to obtain rank-2 and rank-3 solutions. The results of our 
analysis are summarized in Tables 1 and 2. Several of them possess a certain amount 
of freedom. They depend on one or two arbitrary functions of one or two Riemann 
invariants, depending on the case. The range of the types of solutions obtained depends 



on different combinations of the vector fields Xa as given in (1.6). For convenience, we 
denote by EiEj, EiSj, SiSj, EiEjEk, etc, i,j,k = 1,2,3, the solutions which are the 
result of nonlinear superpositions of rank-1 solutions associated with different types of 



wave vectors (1.4 i) and (1.4 ii). By r^,r^ and we denote the Riemann invariants 
which coincide with the group invariants of the differential operators Xa of the solution 
under consideration. 

The arbitrary functions appearing in the solutions listed in Tables 1 and 2 allow 
us to change the geometrical properties of the governed fluid flow in such a way as 
to exclude the presence of singularities. This fact is of special significance since, as 
is well known [161 12] , in most cases, even for arbitrary smooth and sufficiently small 
initial data at t = to the magnitude of the first derivatives of Riemann invariants 
becomes unbounded in some finite time T. Thus, solutions expressible in terms of 
Riemann invariants usually admit a gradient catastrophe. Nevertheless, we have been 
able to show that it is still possible in these cases to construct bounded solutions 
expressed in terms of elliptic functions, through the proper selection of the arbitrary 
functions appearing in the general solution. For this purpose it is useful to select 
DCs corresponding to a certain class of the nonlinear Klein-Gordon equation which is 
known to possess rich families of bounded solutions [1]. We choose elliptic solutions 
of the Klein-Gordon equation because a group theoretical analysis has already been 
performed [22]. The obtained results can be adapted to the isentropic ideal compressible 
fluid flow in (3 -|- 1) dimensions. Thus, we specify the arbitrary function(s) appearing 
in the general solutions listed in Tables 1 and 2, say (f), to the differential constraint 
in the form of the Klein-Gordon 0®-field equation in three independent variables r^, 
and which form the coordinates of the Minkowski space M(l, 2) 



ceM. (2.1) 
Here, we choose to be timelike and r^, to be spacelike coordinates. It is well known 



see e.g. [22]) that equation (2.1) is invariant with respect to the similitude Lie algebra 
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sim(l, 2) involving the following generators 

D = r'drr - ^4>d^, Pi = dr^, i = 1, 2, 3, 

Lab = r^d^t - r^dra, a ^ 6 = 2, 3, (^-^^ 
Kia = -{r^dra - r^dri), a = 2,3, 

where D denotes a dilation, Pi represents translations, Lab stands for rotations and Kia 
for Lorentz boosts. A systematic use of the subgroup structure |22] of the invariance 



group of (2.1) allows us to generate all symmetry variables ^ in terms of the Riemann 
invariants r^,r^,r^. We concentrate here only on the case when symmetry variables 
are invariants of the assumed subgroups Gi of Sim{l,2) having generic orbits of codi- 
mension one. For illustration purposes, we perform a symmetry reduction analysis on 
four selected members of the list of subalgebras given in ([22], Table IV) which involve 
dilations. For each selected subalgebra in Minkowski space M(l,2), we compute the 



group invariants ^ of the corresponding Lie subgroup and reduce the equation (2.1) to 



a second order ODE. The application of the symmetry reduction method to equation 



(2.1) leads to solutions of the form 

0(r) = a(r)F(^(r)), r = {r\ r'^ , r^) , (2.3) 
where the multiplier a{r) and the symmetry variable ^(r) are given explicitly by group 



theoretical considerations and F(^) satisfies an ODE obtained by substituting (2.3) into 



equation (2.1). The results of our computation are listed below. 



1. {D,Pi}: a = {AcUr^Y + (r^f]}-^/\ e = -arctan^, F" + F + = 0, 

2 

2. {D,L,,}: a = {-c(rT/4}"^^ ^ = + 



({1+OF"+(2(+^'\f'+^F + F^ = 0, 



q c 

F" + 1^7^' + = 0, q = -1/3,1- 2,4-3/, / G Z+, 

4. {D + ^-K^2, L, - Kr,} : a = {9/AC)'/\r' - (r^ + r'^/A}-'/', 

^ _ 6{r^ - r^) + 6r\r^ + r^) - (r^ + r^f n ^ p// , I p/ , 1 p , p5 _ n 
^ 8(r2-(ri+r3)V4)3/2 ' U + ^ + 3^ + 3^ + ^ u. 



The parity invariance of (2.1) suggests the substitution 
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which transforms the equations hsted in (2.4) to 
{D,Pi} H" 



i^i2,i:23} H" 

q 



2H 

jj,2 



2{H + H 
1 



2H ^(1 + 



2^+^]H' + ^H + 2H^ 



(2.5) 
(2.6) 



H' 



/2 



2H 



m 
J 



H' + 2H^ 



m 



3q + l 
2q + l 



H' 



/2 



-£H' + -H + 2H^ 
3^ 3 



(0,4/3,2), (2.7) 



(2.^ 



2H 1+^2 

where the three admissible values for the scalar m come from group theoretical consid- 
erations [22]. Each of these four equations possesses a first integral 
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(2.9) 



in which the four sets of functions G,g and constants eo,co obey the respective condi- 
tions 



G 
G 
G 
G 



2 4eo 
-r^ 9 = — , 
4 Co 



3ci 
'T 

3co 
' 4 

3co 
' 4 



64en 



Co 



{a,eo,g') = {0,0,h), {4/3,0,he^'), {2,eo,-^e' 



Co 



3co^^2 
' 4 



(e^ + 1), <? = A:i(l+a'/', eo = 



{ki denotes an arbitrary nonzero real constant). Under a transformation {H,$^) — * {U, () 
which preserves the Painleve property, 

2 

^(0 = uiO/giO, 



G(?2 



the equation (2.9) becomes autonomous 



U'^ - cqU^ - 12eoU^ - AK'U = 0, co ^ 0. 

When K' = 0, [/"^ is either a sine, cosine, hyperbolic sine or a hyperbolic cosine 
function, depending on the signs of the constants, therefore bounded solutions are 
easily characterized. 

When K' ^ 0, it is convenient to first integrate this elliptic equation in terms of the 
Weierstrass function p{(, g2, gs), 

U{0 = -7^ , g2 = I2el = -8el - c,K'\ 

p'2 = 4(p _ ei){p - e2){p - eg) = Ap^ - g2p - ^3, 
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(where we abbreviate p{C, 92,93) by p{C)) then to use the classical formulae which 
connect p and various bounded Jacobi functions. 

This correspondence is quite easy to write down if one uses the symmetric notation of 
Halphen [T3] to represent the Jacobi functions. Halphen introduces three basis functions 

ha{u) = \/piu) - Ca, a = 1,2, 3, 

and the connection between the Weierstrass p function and the Jacobi copolar trio 
cs, ds, ns is given by [14, p. 46] 

cs{z\k) ds{z\k) ns{z\k) u 1 ^2 ^2 — 63 

hi{u) h2{u) h-s{u) z Vei - 63 ' 61-63' 

where k is the modulus of the Jacobi elliptic functions. For full details on Halphen's 
symmetric notation, see [O] . We give here explicit solutions in terms of the Weierstrass 
function, leaving the conversion to Jacobi's notation to the reader. These solutions are 
obtained by convenient choices of the normalization constants cq, cq, K' , k\. 

With the normalization cq = —1/3, cq = —4/3, K' = C, the general solution of 
^ is 

^^K) - ^^^r c - e * ^ 1 93 4 + c e M. (2.10) 



With the normalization cq = k^'^/AS, cq = — (4/3)/cq^, K' = C, the solution of (2.6) 
has the form 



^ , , /- 1 1 4C2 



with fco, C G M. 

The three cases for equation (2.7) associated with the subalgebra {D + ^-^Ki2, L23} 
yield the respective solutions 

q = -k/?,:F\i) = —, ( = ^,32 = 0, 93 



p{i) 3 

g = 4 - 3fc : F\i) = C = Sfco^'/', 92 = 0, 93 = 

q = k — 2 : F^(f) = — — ^ — , C = kr)\os:f, oo = — rr, 9s = H rr- 

PiO-jh I2kl' 216kl 3ki 



2 



Finally, equation (2.8) integrates as (equation no 4 in (2.4)) 



p(C) ^ V2'6'2' V 3A;2 

where 2-^1 denotes the hypergeometric function. 



8 



Using these results, we construct bounded rank-3 solutions of the equations (1.1). 
For this purpose, for each general solution appearing in Tables 2 and 3, we introduce the 
arbitrary functions into the Klein-Gordon equation (2.1) and select only the solutions 
expressed in terms of the Weierstrass p-function. 

For illustration, let us now discuss the case of the rank-3 entropic solution EiE2E^ 
which represents a superposition of three rank-1 entropic solutions Ei given by (1.12). 
We assume that the entropic wave vectors A^^ and A'^^ are linearly independent 
and take the form 

= (a + r ■M,-e^), |e'p = l, i = 1,2,3. 
Hence the corresponding vector fields Xi and Riemann invariants are given by 

d det(e\e^,e^) d Bi d 132 d , , _^ ^, ^ 

where Pi = (e^ x e^)i{a + ■ u) + (e^ x e^)i{a + e ^ ■ m) + (e^ x e^)i{a + ■ u). The 
rank-3 entropic solutions invariant under the vector field X are obtained through the 
change of coordinates 

i = t,x^ = r-^^x, u),x'^ = r'^{x, u),x^ = r^{x, u),a = a,u^ = u^, = u^, = u^, 

(2.12) 

on M.'^ X M^. Specifying the form of the solution as a linear superposition of rank-1 
solutions (1.12), 

a = ai{r^) + a2{r^) + as{r^), u = K{e^ai{r^) + e^ttiir^) + e^a^ir^)) (2.13) 



the fluid dynamics equations (1.1) transform to 

2 3 

J2 E ['^(^"^ ■ ^"'')' + (1 - '^)(^"'^ ■ ^"'') - ^Kir'mr') = 0, (2.14) 

1=1 j=i+l 

while the invariance conditions have the form 

flj- = u} = uf = ul = 0. (2.15) 

This solution exists if and only if the three entropic wave vectors e^, e^, intersect at 
a certain specific angle given by 

cos(/)ij = --, ii7^j = l,2,3, 

where (pij denotes the angle between the wave vectors and flUl [20] . Imposing the 
condition that each of the functions cij(r*), z = 1, 2, 3 obeys the ODE 



F" + F + F^ = 



9 



then according to (2.10), the rank-3 entropic solution takes the form 



i=i [pij"", |, ^ + \Cf) + I) ' i=i [p{r\ |, ^ + \Cf) + I) ^ 

-(l + /«) —t + X-x, z = 1,2,3. 

(p(rS|,A + |C^) + l)V^ 



(2.16) 



Making use of an exphcit expression for the zeros of the p-function, we show that the 



values for which the denominator in solution (2.16) vanish are not located on the real 
axis for a specific choice of the constants of integration Ci. Then we have the following 
result. 



If the constants of integration Cj are equal to v 19/6, then the elliptic rank-3 entropic 



solution (2.16) of the isentropic ideal compressible fluid flow equations (1.1) is bounded. 

Indeed, according to recent results obtained by Duke and Imamoglu in j6|, the 
location of the zeros of the p-function can be given explicitly in terms of generalized 
hypergeometric functions. 

Considering a lattice £ = Z + rZ, Imr > 0, the doubly periodic Weierstrass p- 
function is defined by 



where the sum ranges over all u E C. Note that the p-function assumes every value of 
the extended complex plane exactly twice in C and since it is even, its zeros are of the 
form ±^0- The value of zq can be determined from the following theorem 0. 

The zeros ±zq of the p-function are given by 

where pFq denotes the generalized hypergeometric functions and s is a function of the 
modular discriminant A and the Eisenstein series E^^, 



s=l-'-^, Air) = qllil - q-r\ i54(r) = 1 + 240 q = e'^ . 

(2.18) 

It is understood that the principal branch is to be taken in the radical expression s^^^. 

We now illustrate the application of this theorem by showing that the function 
pjz, 4 /3, 1) is always strictly positive for real z. This case is the specific case of solution 
(2.10) for which C = VT9/Q. From given invariants g2 and gs, by finding the roots 
61,62,63 of the cubic polynomial 4t^ — g2t — g^, one can determine the values of the 
periods ui and 002- In the case when g2 = 4/3 and g^ = 1, we obtain the periods 
iui = 2.81 and UJ2 = 1.405 + 12.902, hence r = c^a/c^i = 0.5 + il.033. 
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For any pair of periods 001,002, the lattice wiZ + L02IJ can be rescaled in such a way 
that LOi is normahzed to 1 by using the well-known formula for the p-function 



p{z,UJi,UJ2) = p{z/uJi,UJ2/uJl)/ujl. 



(2.19) 



Introducing the numeric value of r into (2.18), we can evaluate from (2.17) and (2.19) 
the zeros of the Weierstrass function p{z, 4/3,1), zq = ±(1.405 + iO.929), which are 
indeed complex. Note that the p-function always possesses a double pole at 2; = and 
that it tends to +00 at this point. Since it is continuous for all real z G {0,ooi), this 
implies that it must always be strictly positive. Hence, p{z, 4/3, 1) + 1/3 > on the real 
interval (0,1^71) and the function (p(z,4/3, l) + l/3)~^ is therefore bounded for all real z. 



Choosing the constants Cj as values of the initial constant C for which solution (2.10) 



is bounded then guarantees the boundedness of functions a and u in solution (2.16) 



Therefore, solution (2.16) is bounded everywhere, even when the Riemann invariants 
r* admit the gradient catastrophe. A similar analysis can be applied to every solution 
presented in Table 3 to show that they are bounded. This can be accomplished by the 
same procedure as presented above through the use of the theorem from [6]. 

This solution is physically interesting since it remains bounded for every value of 
the Riemann invariants r*. Thus, it represents a bounded solution with periodic flow 
velocities (see figure [T|. Similarly, it is possible to submit the arbitrary functions of the 
differential equations no 2, 3, 4 listed in (2.4 ) to obtain other types of bounded solutions. 
Table 3 presents these various types of solutions with their corresponding Riemann 
invariants. They are all bounded solutions of periodic, bump or kink type. Note that 



these solutions of (1.1) admit gradient catastrophes at some finite time. Hence, some 



discontinuities can occur like shock waves [5l |T3]. Note also that the solutions remain 
bounded even when the first derivatives of r' tend to infinity after some finite time T. 
However, after time T, the solution cannot be represented in parametric form by the 
Riemann invariants and ceases to exist. 



3 Concluding remarks 

The methods presented in this paper can be applied quite broadly and can usually 
provide at least certain particular solutions of hydrodynamic type equations. The 
conditional symmetries refer to the symmetries of the overdetermined system obtained 



by subjecting the original system (1.1) to certain differential constraints defined by 



setting the characteristics of the vector fields Xa to zero. The conditional symmetries 



are not symmetries of the original system (1.1). However, they are used to construct 



classes of rank-3 solutions of this system which are not obtainable by the classical 
symmetry approach. Among the new results obtained, we have rank-2 and rank-3 
periodic solutions expressed in terms of the Weierstrass p-function that we have shown 
to be bounded over the real axis. They represent bumps, kinks and multiple-wave 
solutions, all of which depend on Riemann invariants. These solutions remain bounded 
even when the invariants admit the gradient catastrophe. 
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Figure 1: Bounded entropic rank-3 solution 



Among the questions that one may ask is what role do exact analytical solutions 
play in the physical interpretation. One possible response is that such solutions may 
display qualitative behaviour which would otherwise be difficult to detect numerically 
or by approximations. For example, the doubly periodic properties of certain solu- 
tions expressed in terms of the Weierstrass p-function would not be very easily seen 
numerically. 

One could also inquire about the stability property of the obtained solutions. In- 
deed, solutions which possess the property of stability should be observable physically 
and such analysis could be the starting point for perturbative computations. This task 
will be undertaken in a future work. 
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